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Method of analyzing a data set comprising a volumetric representation of an object to be 
examined 



The present invention relates to a method of analyzing a data set of an object 
to be examined, which data set comprises voxels of at least a first type and a second type. 

Such a volumetric representation of the object can be acquired by means of 
various techniques. Some examples thereof are: 3D X-ray rotational angiography, computed 
5 tomography, magnetic resonance imaging or magnetic resonance angiography. The object 
typically represents a patient to be examined. 

The international patent application EPOO/09505 (PHN 17.678) of the same 
applicant (not yet published) relates to a method of the type mentioned above. This method 
concerns the determination of the volume of an aneurysm in a blood vessel. According to the 
10 method described in the cited international patent application EPOO/09505, a self-adjusting 
probe is defined for analysis of the data set. Doing so requires a lot of expert interaction, and 
is a burden since the analysis is often performed real-time during vascular surgery. 

It is an object of the invention to provide a method of the type described above 
1 5 which requires less expert interaction. 

To this end, the method according to the invention comprises the following 

steps: 

a) classifying the voxels as voxels of the first, the second or further types; 
20 b) determining which of the voxels of the first type are boundary voxels that adjoin voxels 
of the second or further types; 

c) assigning a data value to each voxel of the first type, said data value representing a 
measure of the distance between said voxel and the nearest boundary voxel; 

d) classifying the voxels of the first type that have a distance data value exceeding a 
25 predetermined threshold as aberration voxels. 

The method according to the invention allows for quick semi-automatic 
determination of the voxels defining the aberration. 

A preferred version of the method also comprises the following steps: 
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e) determining which of the aberration voxels are boundary aberration voxels adjoining 
voxels of the first type; 

f) adding a number of voxels of the first type that form a shell of a certain thickness to the 
aberration voxels. 

5 The number of voxels defining the aberration can now be determined with a 

greater accuracy. 

Preferably the step f comprises the steps of: 
fl . assigning a data value to each voxel of the first type, said data value representing a 
measure of the distance between said voxel and the nearest boundary aberration voxel; 
1 0 f2. classifying the voxels of the first type that have a distance data value less than or equal to 
a predetermined ceiling value as aberration voxels. 

In one version the volume of the aberration is computed by determining the 
sum of all aberration voxels and multiplying the sum of the aberration voxels by the volume 
of a single voxel. 

15 In an elegant version of the method according to the invention said distance 

data values are computed by means of a distance transform function. 

It is to be noted that the distance transform function is knovra per se in a field 
other than that of the present invention, that is the field of computer vision where the distance 
transform function is used for object recognition based on pictures. The distance transform 
20 function is described, for example, in the paper "Distance transformations in arbitrary 

dimensions" by G. Borgefors, Computer Vision, Graphics and Image Processing, Vol. 27, 
pages 321-345, 1984. 

The threshold and/or ceiling value can be set by the user. In order to reduce 
the necessary expert interaction even further, the threshold and/or ceiling value can 
25 alternatively be computed on the basis of a histogram of distance data values. 

To further enhance the accuracy of the computation of the voliraie of the 
aberration, another version of the method according to the invention comprises the steps of: 

defining a tubular structure of voxels of the furst type piercing through the 
aberration; determining the nvimber of voxels of the tubular structure and subtracting said 
30 number from the number of aberration voxels. Preferably, the method also comprises the 
steps of: 

a. classifying all boundary aberration voxels as potential tubular structure voxels; 

b. selecting a starting point from among the potential tubular structure voxels; 

c. selecting an end point from among the potential tubular structure voxels; and 
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d. connecting the starting point to the end point, thus defining the tubular structure. 

The invention also relates to a computer program for carrying out the method 
according to the invention. 

5 

The invention will be described in detail hereinafter with reference to the 
attached drawings, in which: 

Figure 1 shows schematically a volumetric representation of an object to be 
examined after application of a first version of the method according tot the invention, and 
10 Figure 2 shows the volumetric representation of figure 1 after application of a 

second version of the method according tot the invention. 

Figure 1 shows an example of a possible result of application of the present 
invention to an object to be examined. In this case the object is part of the vascular system of 
15 a patient. Using of well-known techniques the object is recorded and represented by means of 
a data set comprising volume elements (voxels). An example of such a technique is described 
, for example, in the article "3D rotational angiography: Clinical value in endovascular 
treatment", by Moret et al., published in Medicamundi, Vol. 42, Issue 3, November 1998. In 
the present example there are generally two kinds of voxels, that is, voxels of the fnst type, 
20 which in this case are vessel voxels, and voxels of the second type, which in this case are 
tissue voxels (not shown). 

Figure 1 shows schematically the object 1 which comprises blood vessels 2 
that connect with an aneurysm 3, being is an aberration of the vessel acting as a reservoir. 

Application of the method according to the invention allows for determination 
25 of the volume of the aneurysm, this is essential for successful and safe medical treatment 
thereof. 

In a first step vessel voxels are classified as voxels of the first type and tissue 
voxels are classified as voxels of the second type. This operation can be performed by means 
of known techniques, such as a 'region growing algorithm'. 
30 In a next step it is determined which of the vessel voxels adjoin tissue voxels. 

These voxels are referred to as boimdary vessel voxels. Preferably, this determination is 
performed by finding the facing neighbors, that is neighbors which have one voxel face in 
common. 
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Subsequently a data value representing a measure of the distance between the 
relevant vessel voxel and the nearest boundary vessel voxel is assigned to each vessel voxel. 
According to the invention said distance data values are computed by means of a distance 
transform function. 

5 According to this distance transform function the distance data values of the 

boundary vessel voxels are set at a low value, for example, e.g. zero. For all other voxels 
initially a high value is chosen, for example 127. The distance transform function then 
computes distance data values for each voxel. An elegant way of computing the distance data 
value is as follows: the distance data value for a specific voxel is defined as the minimum of 

1 0 the distance data value of each of the neighbors of that specific voxel plus a number, for 

example, one. Preferably the distance transform fimction is computed once for all voxels in a 
certain order. Next the distance transform is again computed for all voxels, but now in the 
reverse order. In order to simplify and hence shorten the computation preferably only vessel 
voxels are involved in the computation. An example of the result of this computation is 

15 shown in figure 1 . 

All vessel voxels with a distance data value exceeding a predetermined 
threshold are then classified as aberration voxels. The user can set this threshold 
interactively. Preferably the threshold should be chosen to be slightly greater than the radii of 
the connected vessels 2. As an altemative the threshold can be computed automatically, for 

20 example, on the basis of a histogram of distance data values. In this example let it be 

assumed that the threshold is three. All vessel voxels with a distance data value equal to or 
larger than three are then classified as aberration voxels. In the present example the 
aberration voxels are aneurysm voxels and a first estimate for the number of aneurysm voxels 
is now available. 

25 In order to compute the number of aneurysm voxels more precisely the next 

steps are performed. 

These steps mainly mean that a new distance data value is computed for a 
select number of voxels to make sure that all boundary aberration voxels are also included. 
Again the distance transform fimction can be used for this purpose. For useful techniques and 
30 remarks on performing these steps we refer to those mentioned above. 

First it is determined which of the aberration voxels are boundary aberration 
voxels adjoining vessel voxels. 
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Next a data value representing a measure of the distance between the reference 
vessel voxel and the nearest boundary aberration voxel is assigned to each vessel voxel. The 
result of this computation is shown in figure 2. 

All vessel voxels with a distance data value less than or equal to a 
predetermined ceiling value are now classified as aberration voxels. Preferably, the ceiling 
value is equal to or slightly greater (e.g. by one unit) than the threshold. In this case the 
ceiling value is equal to the threshold, that is, three. As can be seen when comparing figure 1 
and figure 2 the second computation results in a more accurate determination of the number 
of aneurysm voxels. 

Now that the aneurysm voxels are known they can simply be summed and 
multiplied by the volume of a single voxel so as to determine the volume of the aneurysm 3. 

In the specific case of determining the volume of an aneurysm those voxels 
that represent part of the vessel which is needed for an unrestricted blood flow between the 
vessels 2 should preferably not be included. To account for these voxels a tubular structure of 
vessel voxels piercing the aneurysm is defined. 

This can be performed by means of a method comprising the following steps. 
First all boundary aneurysm voxels are classified as potential tubular structure voxels. Out of 
the potential tubular structure voxels one voxel is selected as a starting point and one voxel is 
selected as an end point. To define the tubular structure the starting point and the end point 
are connected. This connection can be performed by means of known techniques, such as the 
so-called "self adjusting probes" mentioned above. 

On the basis of the definition of the tubular structure the corresponding 
number of voxels can be determined and subtracted from the number of aneurysm voxels. 

On the basis of the above explanation of the method of the invention a skilled 
person will be able to translate the steps of the method into a computer program for carrying 
out the method. 

Summarizing the invention refers to a method for analyzing a data set 
comprising a volumetric representation of (a part of) an object to be examined in order to 
detect a sub-object of a volume greater than the volume of surrounding sub objects and to 
measure the volume of said sub object. This method is especially useful in the field of 
medical diagnostics and treatment where the object notably is a patient to be examined. 

The invention is, of course not limited to the versions described or shown, but 
generally extends to any version that is, which falls within the scope of the appended claims 
as seen in the light of the foregoing description and drawings 



